Three-Loop Radiative-Recoil Corrections to Hyperfine Splitting 
Generated by One-Loop Fermion Factors 



Michael I. Eides * 

Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA 
and Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188300, Russia 

Howard Grotch''' 

Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA 

Valery A. Shelyuto ^ 

D. L Mendeleev Institute of Metrology, St. Petersburg 198005, Russia 



Abstract 

We consider three-loop radiative-recoil corrections to hyperfine spUtting in 
muonium generated by diagrams with one-loop radiative photon insertions 
both in the electron and muon lines. An analytic result for these nonloga- 
rithmic corrections of order a{Z'^a){Za){m/M)Ep is obtained. This result 
constitutes a next step in the implementation of the program of reduction of 
the theoretical uncertainty of hyperfine splitting below 10 Hz. 

I. INTRODUCTION 

Three-loop radiative-recoil corrections to hyperfine splitting in muonium are enhanced 
by the large logarithm cubed of the electron-muon mass ratio [1] (see, also review [2]). 
The leading logarithm cubed contribution is generated by the graphs with insertions of the 
electron one-loop polarization operators in the two-photon exchange graphs. It may be 
obtained almost without any calculations by substituting the effective charge a{M) in the 
leading recoil correction of order [Z a){m / M)E p and expanding the resulting expression in 
a power series in a. Calculation of the logarithm squared term of order a'^{Za){rn/M)Ep 
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is more challenging [3]. Different graphs generate logarithm squared terms, and all such 
contributions were obtained a long time ago [1,3,4]. The sum of the logarithm cubed and 
logarithm squared terms is given by the expression ^ 

A£;=f-lln3^ + lln^^)^^E.. (2) 

Due to recent experimental and theoretical progress, single-logarithmic and nonloga- 
rithmic contributions of orders a'^{Za){m/M)EF and a{Z'^a){Za){m/M)EF to hyperfine 
splitting in muonium are now also phenomenologically relevant. Numerous sets of gauge 
invariant diagrams generate single-logarithmic and nonlogarithmic contributions. 



fi 1-loop \ fi 1-loop I 

^•\l-loop S 1-loop S 



FIG. 1. Diagrams with two fermion factors 

Below we consider three-loop radiative-recoil corrections to hyperfine splitting in muo- 
nium generated by the diagrams in Fig. 1. These diagrams are obtained from the skeleton 
diagrams in Fig. 2 by making all one-loop radiative photon insertions both in the electron 
and muon lines. The two-loop radiative-recoil corrections generated by all one-loop radiative 
photon insertions only in the electron or only in the muon line were calculated earlier (see, 
e. g., review in [2]). The leading recoil correction of order Za{m/ M)Ep is generated by the 
skeleton diagrams with two exchanged photons in Fig. 2. The characteristic loop momenta 
in the skeleton diagrams are larger than the electron mass, and therefore the leading recoil 
correction to hyperfine splitting may be calculated in the scattering approximation, ignoring 
the wave function momenta of order mZa (see, e. g., [2]). It was obtained a long time ago 
[5-7] 



3mM Za , M ~ 

^^rec = J— In — Ef. 3 



We define the Fermi energy as 

where m and M are the electron and muon masses, rrir is the reduced mass of the electron- muon 
system, a is the fine structure constant, c is the velocity of light, h is the Planck constant, Roo 
is the Rydberg constant, and Z is the nucleus charge in terms of the electron charge {Z = 1 for 
muonium). The Fermi energy Ep does not include the muon anomalous magnetic moment 
which does not factorize in the case of recoil corrections, and should be considered on the same 
grounds as other corrections to hyperfine splitting. 
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Note that this correction is proportional to the logarithm of the electron-muon mass ratio, 
and quite remarkably it turns out that the logarithmic contribution is a complete result, 
there is no nonlogarithmic contribution of this order. 

The radiative insertions can only increase the characteristic integration momenta in the 
diagrams in Fig. 2 and hence the scattering approximation remains valid for calculation of 
two- and three-loop radiativc-rccoil corrections. The two-loop radiative-recoil corrections of 
order a{Z a){m / M)E p generated by the one-loop radiative photon insertions in the electron 
line are logarithmic in the electron-muon mass ratio. Since the leading recoil correction 
of order Za{m/M)EF is hnear in the logarithm of the large mass ratio, one could expect 
that the correction of order a{Za){m/M)Ep is proportional to the logarithm squared. This 
does not happen and the logarithm squared contributions cancel as was first discovered 
in [8] by direct calculation. The simplest way to understand this cancellation is to recall 
that in the Landau gauge radiative insertions in the electron line are nonlogarithmic [9], 
and hence, being gauge invariant, the sum of these insertions is nonlogarithmic in any 
gauge. The single-logarithmic term of order a{Za){m/M)Ep was obtained in [10], and the 
nonlogarithmic terms were calculated numerically in [11] and analytically in [12] 
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(4) 



One more feature of the calculations in [11,12] deserves to be mentioned. One-loop radiative 
insertions in the electron line include the terms connected with the one-loop anomalous 
magnetic moment. These terms have different low energy behavior in comparison with all 
other terms in the dressed electron line and could in principle complicate calculation of the 
radiative-recoil corrections. However, as was discovered in [11,12] the terms connected with 
the one-loop anomalous magnetic moment do not give any contribution at all to the radiative- 
recoil corrections of order a{Za){m/M)Ep. Finally, let us mention that numerically the 
nonlogarithmic part of the correction of order a{Za){m/M)Ep is rather large, of order tt^, 
which is just what one should expect for the constants accompanying the large logarithm. 

The two-loop radiative-recoil corrections of order {Z'^a){Za){m/M)Ep are generated by 
all one-loop radiative photon insertions only in the muon line in the diagrams in Fig. 2, and 
were obtained in [11,13] 



Q 3Q 
^C(3)-37rMn2 + | 



[Z^a){Za)m ~ 

Hip. 
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(5) 



Two features of these corrections deserve to be mentioned. First, radiative insertions in the 
muon line do not generate logarithmic terms at all, as can be understood with the help of 
the generalized low energy theorem [10,14]. Second, just as in the case of the insertions in 
the electron line the terms connected with the muon anomalous magnetic moment do not 
give any contribution to the radiative-recoil corrections of order {Z'^a){ZQ)(rri/M)Ep. 

In this work we analytically calculate three-loop radiative-recoil corrections to hyperfine 
splitting in muonium generated by the diagrams in Fig. 1 with all one-loop radiative photon 
insertions both in the electron and muon lines. We show that these corrections are nonlog- 
arithmic and unlike the case of the radiative-recoil corrections of orders a{Za){m/M)Ep 
and {Z^a){Za){m/M)Ep the one-loop anomalous magnetic moments of both particles give 
nonvanishing contributions to the correction under investigation. 
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FIG. 2. Diagrams with two-photon exchanges 

II. GAUGE INVARIANT REPRESENTATION FOR RADIATIVE 

CORRECTIONS 




Let us consider the general structure of the radiative-recoil corrections in Fig. 1. To 
this end it is convenient to introduce the one-loop fermion factor L^i,{k) as a sum of the 
diagrams in Fig. 3. In terms of the electron and muon factors the radiative- recoil contribution 
to hyperfine splitting generated by the ladder and crossed-ladder diagrams in Fig. 1 has the 
form 



TT 



(6) 



The sum of the electron factors CjfJ{k) = L^^},{k) + A;) which enters eq.(6) for the 

radiative corrections is just the gauge invariant Compton scattering amplitude for a vir- 
tual photon, and satisfies the identity k'^Cj^J{k) = 0. The electron Compton amplitude 
is invariant under the substitution k ^ —k and /i ^ u, and hence, we can substi- 
tute the muon Compton amplitude instead of the muon factor in the integral in eq.(6) 
L^^J{—k) — >• [Ljf^)(— /c) -I- Ll^^{k)]/2 = C^^\—k)/2 obtaining a more symmetric expression 
for the energy shift 
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FIG. 3. Fermion factor 

To simplify further calculations we represent the electron and muon Compton amplitudes 
in eq.(7) as sums of two gauge invariant terms (we write the formula only for the electron, and 
the respective expression for the muon is obtained by the substitution m — > M, a — > Z'^a) 



where 



C]^f\k) 



a 1 
27r2m 



k"^ — 2'mko 



(8) 



(9) 



crup{-k'')(p + k + m)j^ + j^(p + k + mjai^pk'' 
k"^ -\- 2mko 
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corresponds to the anomalous magnetic moment, and C^^jf>{k) includes all other terms. 

It is easy to check directly that Cj^j"'\k), and hence, C^^J'\k) are gauge invariant. The 
breakdown in eq.(8) is helpful because Cj^j'^^k), and Cjfj^\k) have different behavior at small 
photon momenta k. As we will see below this different low energy behavior determines the 
structure of integrals for the contributions to hyperfine splitting. 

We can further simplify the amplitude Cjfj^^k), preserving only the terms which con- 
tribute to hyperfine splitting. The simplified expression (still satisfying the Ward identity 
k^'Cj^f^k) ^ 0) has the form 

^^^^"'-^ k^-Aw?kl ■ 

In terms of the representation in eq.(8) the contribution to hyperfine splitting in eq.(6) 
can be written as a sum of three gauge invariant terms 



. ^ 3 {Za)mM ~ f d^k 
/\E — Ef 



16 TT J m'^k^ 



Cif\k)C\l^f\-k) + &\k)C^^f\-k) (11) 



+C(^j(A;)C(^'^)(-A;)] = /^E^ + /^E^^ + /^E^^^ . 

It is important to note that we know in advance that there is no logarithm of the mass 
ratio in the sum of all contributions in eq.(ll). Such a logarithm can only arise from the 
integration region m < k < M, where the electron factor is in the asymptotic regime. The 
asymptotic expression for the electron factor was calculated e.g., in [3], and contains only 
the skeleton spinor structure j^kju. On the other hand, all terms in the muon factor except 
the term with the muon anomalous magnetic moment are, in this integration region, addi- 
tionally suppressed by an extra factor /c^/M^ in comparison with the logarithmic skeleton 
integral, and thus cannot produce a logarithmic contribution. As to the term with the muon 
anomalous magnetic moment, its contribution to the recoil integral vanishes identically due 
to its spinor structure, see, e. g., [14]. 

III. TWO ANOMALOUS MAGNETIC MOMENTS 

Let us start our calculation with the term AE^ connected with the product of two one- 
loop anomalous magnetic moments. Projecting the spinor structures of the fermion factors, 
written in the form of eq.(9) or eq.(lO), on hyperfine splitting we obtain 

^p/_ ^(^^n)(Zn)mM ~ r d'k k^k ^ + ik^ + hj) 

TT^ in^k^ {k"" - Aklw?){k'' - AklM^y ^ ' 

or after the Wick rotation and transition to the spherical coordinates /cq = kcos9, |k| = 
ksinO 

TT^ Jo Jo (A;2 + 4m2 cos2 9) {k^ + AM^ cos^ 9) ^ ' 
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Calculating the angular integral wc discover that the remaining momentum integral diverges 
like dk'^/k^. This divergence indicates the existence of the nonrecoil correction of order 
a{Z'^a), which is of lower order in Za. It is connected with the one-photon exchange, 
and is well known. We subtract this power divergence and, after the subtraction, obtain a 
convergent integral, which can be easily calculated 



^ 9 a{Z^a){Za) mM ^^M ~ 



IV. SUBTRACTED ELECTRON FACTOR AND THE MUON ANOMALOUS 

MAGNETIC MOMENT 

The second contribution in eq.(ll) arises from the product of the subtracted electron 
factor and the muon anomalous magnetic moment, and we write it in the form 

^^j, ^ 3 iZa)mM ~ f d'k 
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f^&^j (15) 

The muon Compton amplitude is gauge invariant and satisfies the Ward identity 
k^C^Ji;"'\k) ~ 0. Therefore, we can omit all terms in the subtracted electron factor which are 
proportional to A;^. This means that we can use the expression for the subtracted electron 
factor from [15,16], where all terms proportional to k^ are thrown away. We represent this 
electron factor as a sum of seven terms L^^J^^{k) — Y^lL^/llik), which are 

L^m + - /; d^i' + c..^) . (16) 



k'^ + 2mbko + m?a'^ \ m 



where A;^ = /cq — k^. Each term in the electron factor corresponds to the respective coefficient 
function q , and explicit expressions for these coefficient functions are collected in Table I. 

We preserve in the electron factor only the spinor structures (7^fc7jy)(e) and {j^ju){e) relevant 
for hyperfine splitting, and the projection on hyperfine splitting is understood. Auxiliary 
functions of the Feynman parameters a{x, y) and b{x, y) are defined by the relationships 

a^^ 6=^. (17) 
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TABLE I. CoefHcients in the Fermion Factor 



Cl 


16 


{l — x){x — Zy) — 2ylna; 


!/(l-2/)^ 


C2 


4 


- (1 - x) (x - y - 2y2/x) + 2{x - 4y + Ay'^/x) In x 


2/(1-2/)'' 


C3 


1 


1 - 6x - - (y/x)(26 - 6?//x - 37x - 2^^ + 12x|y + 161nx) 


2/(1-2/)^ 


C4 




2x - - 5y + 7xy ^ 


2/(1-2/)^ 


C5 




6x — 2>x^ — 8y + 2xy j 


2/(1-2/)^ 


C6 


a;^(l-2/)'' 


C7 


X 
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The explicit expression for the muon factor can be obtained from the expression for the 
electron factor by the substitutions m M and a — > Z'^a. 

Taking projection on the hyperfine splitting and contracting the Lorentz indices, we 
obtain the integral for the contribution to the hyperfine splitting as 



1 

— - / dx dy 

STT-^/X Jo Jo 



d^k 



1 



(18) 



{( 



Qkl 



cik^ + C2{k 



2\2 



-k^ + 2bko + a2)3 



+ 



c^k'^ + C42kQ 
-A;2 + 2bko + a2)2 



-(6 + 4^^0 



c^k'^ + C6p2ko 



-k^ + 2bko + a^y -k^ + 2bko + a^ 



where ji = m/ (2M ) and we rescaled the integration momentum, so that now it is measured 
in units of the electron mass. 

The analytic calculation of the integrals in eq.(18) is one of the more tedious steps in the 
present paper. These integrals are of the same type as the integrals in [15,16], and we use 
for calculations the same methods as in those papers. First we represent each integral as a 
sum of /^-dependent and /x-independent integrals. The ^-independent integrals admit direct 
analytic calculation. To calculate the /x-dcpcndcnt integrals we separate the contributions 
of large and small integration momenta with the help of an auxiliary parameter a such 
that 1 <S (7 <S l//i. In the region of small momenta we use the condition /xfc <^ 1 to 
simplify the integrand, and in the region of large momenta the same goal is achieved with 
the help of the condition k ^ 1. Finally, for k ^ a both conditions on the integration 
momenta hold simultaneously, so in the sum of low-momenta and high-momenta integrals 
all cr-dependent terms cancel, and we obtain a a-independent result for the integral (for 
more detailed exposition of this method see, e.g., [17]). Here we skip the calculations and 
present only the final result 
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a{Z'^a){Za) m 



M 



E 



F- 



(19) 



V. TOTAL ELECTRON FACTOR AND THE SUBTRACTED MUON FACTOR 

Consider now the last contribution 

A^'" = J ^C<S(.)Cr (-^. (20) 

Due to the generalized low energy theorem the subtracted virtual muon Compton amplitude 
Cjj^''^\—k) is suppressed like k'^/M'^ for momenta k < M (see, e.g., [14]). Hence, the recoil 
correction of first order in the small mass ratio arises from the integration region in eq.(20) 
where characteristic momenta are of order M. At these high integration momenta only the 
leading term in the ultraviolet asymptotic expansion of the one-loop electron factor survives 
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in the integral. This leading term was calculated in [3], and up to the terms proportional to 
k/j, and/or ki, has the form 

|2|Z.. (21) 

Due to gauge invariance of the subtracted muon factor k^C^J^^j^^—k) = k^Cjj^'^\—k) = 0, 
and then the terms in asymptotic expansion of the electron factor which are linear in /c^ 
and/or k„ do not give a contribution to the energy shift in eq.(20). Only the term in 
eq.(21) is relevant for the calculation of the leading recoil correction. In addition further 
simplifications can be made. The subtracted muon Compton amplitude also can be written 
as a sum of terms linear in A;^ and/or k^, and the remaining terms. But it is easy to see 
that k^j^kji, — k'^^fi has zero projection on hyperfine splitting, and hence we can omit all 
terms proportional to k^ and/or k^, in the expression for the subtracted muon Compton 
amplitude in eq.(6). Then the radiative- recoil contribution to hyperfine splitting of order 
a{Z'^a){Za)EF in eq.(6) reduces to 

This last integral is proportional to the integral for radiative-recoil corrections of order 
{Z^a){Za){'m/M)Ep generated by radiative insertions in the muon line [13]. Let us recall 
that the leading term in the asymptotic expansion of the skeleton virtual Compton amplitude 
is 

C'(e,.feeO(^) ^ -2^^- (23) 

Comparing this asymptotics with the expression in eq.(21) and using the result of [13] we 
obtain 



AE 



III 



45,, , 157r2 195 

C 3 + n2 

8 ^ ^ ^ 4 32 



a{Z'^a){Za) m 
7r3 M 



Ep. (24) 



VI. SUMMARY 



The total three-loop radiative-recoil correction to hyperfine splitting in muonium gener- 
ated by the diagrams in Fig. 1 with one-loop radiative photon insertions both in the electron 
and muon lines is given by the sum of the contributions in eq.(14), eq.(19), and eq.(24) 



AEi 



15,,^, 157r\ ^ 277r2 147 

— C(3) H In 2+ 

8 ^ 4 16 32 



a{Z'^a){Za) m 



M 



Ep 



(25) 



Note that, as explained in Section II, single-logarithmic contributions cancelled in this result. 
We also would like to emphasize that unlike the case of the radiative-recoil corrections 
generated by the radiative photon insertions only in the electron or only in the muon line, 
the one-loop anomalous magnetic moments of both particles give nonvanishing contributions 
in eq.(25). 
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Some other three-loop radiative-recoil single-logarithmic and nonlogarithmic radiative- 
recoil corrections were also calculated recently. The corrections of order a^{Za){m/M)Ep 
generated by the graphs with two-loop polarization insertions (irreducible and reducible) in 
the two- photon exchange diagrams were obtained in [18] 



AE-: 



6C(3) + 



33 



In 



M 97 



m 



C(3) - 16Li4 



137r2 



8 ' V2 

44951 a^(Za) m 



+ 



In' 2 



ln^2 



E 



36 36 432 J M ^^^^ 

Single-logarithmic and nonlogarithmic corrections generated by the diagrams with one- 
loop polarization insertions in the exchanged photons and radiative photon insertions in the 
fermion lines were obtained in [19]. These are corrections of orders a^{Za){m/ M)Ef and 
a{Z'^a){Za){m/M)EF, and they have the form 



AE. 



22 M 

-In— + 7.36110 (3) 
3 m 



a^{Za) m 



M 



E, 



(27) 



+ 



6C(3) 



47r2 In 2 + — In — + 22.51939(5) 

2 / TTl 



a{Z'^a){Za) m ~ 
Ep, 



TT^ M 

Combining all three-loop single-logarithmic and nonlogarithmic corrections to hyperfine 
splitting in eq.(25), eq.(26), and eq.(27) we obtain {Z — 1 below) 

67\ , M /1\ 2%^ 

12. 



AE; 



tot 



-Att^ In 2 



In— -14C(3)-16Li4fJ 
m \2 



+ 



111^2 



(28) 



157r2, 
H — ln2 



2,4 5n^ 19l7r2 

- In^ 2 + + 

3 36 144 



12959 
864 



+ 29.88049 (6) 



m 



TT 



3M 



Ef 



or 



67 



M 



47r' In 2 + — In — + 45.0546 



12/ m 

Numerically this contribution to the muonium HFS is 

AEtot = -0.019 2 kHz. 



7r3 M 



(29) 



(30) 



Currently the theoretical accuracy of hyperfine splitting in muonium is about 70 Hz. A 
realistic goal is to reduce this uncertainty below 10 Hz (see a more detailed discussion in 
[2,19]). The new contribution obtained in this paper, eq.(25), together with the results of 
other recent research [16,18-22] makes achievement of this goal closer. Phenomenologically, 
the improved accuracy of the theory of hyperfine splitting would lead to a reduction of the 
uncertainty of the value of the electron-muon mass ratio derived from the experimental data 
[23] on hyperfine splitting (see, e.g., reviews in [2,24]). 
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